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INTRODUCTION 


Maximum  overpressures  developed  in  internal  explosions  have  been  the  object  of  study 
by  the  Naval  Postgraduate  School  (NPS)  for  a  number  of  years.  Previous  reports  have 
appeared  on  the  effects  with  a  variety  of  C-H-N-O  fuels  in  air1  and  for  magnesium  alone,2  as 
well  as  magnesium-fuel  mixtures  in  air.3 

The  present  study  represents  the  extension  of  this  effort  to  the  adiabatic  properties 
of  the  internal  explosions  of  aluminized  fuels  in  air. 

Computations  have  been  carried  out  using  the  HP  9845  computer.  The  expanded 
capability,  relative  to  the  WANG  720C  used  in  the  earlier  work,1'3  has  permitted  a  much 
more  detailed  and  precise  analysis  of  this  system,  in  spite  of  its  considerably  greater 
complexity. 

Aluminum  and  fuel  are  considered  to  be  introduced  at  25°C  into  the  air  volume 
at  1  atmosphere.  It  is  assumed  that  combustion  is  so  rapid  that  the  process  can  be  treated 
as  adiabatic  and  that  the  products  are  distributed  uniformly  throughout  the  volume,  which 
is  considered  constant.  The  volume  of  condensed  phases  is  neglected. 

The  fuels  considered  are  listed  in  Table  1 ,  along  with  their  chemical  formulas  and 
a  few  pertinent  properties.  The  “relative  oxygen  content”  given  in  the  last  column  is  a  rough 
means  of  comparing  the  fuels  with  one  another.  It  is  used  as  the  ordinate  in  several  of  the 
figures. 

In  all  the  following,  the  symbol  C  will  represent  the  charge,  or  the  mass  of  fuel,  per 
cubic  meter  of  air  and  M,  the  mass  of  metal  (aluminum)  per  cubic  meter  of  air.  Charge-to- 
metal  (C/M)  ratios  and  total  concentrations  (C  +  M,  in  kilograms  per  cubic  meter)  were 
varied  over  the  range  0.1  to  10  for  each. 

In  several  of  the  tables  a  code  has  been  used  to  indicate  the  condensed  phases  present 
in  the  final  system.  The  key  to  this  code  is  given  in  Table  2.  Again,  in  most  of  the  figures, 
the  curves  have  been  dotted  in  those  regions  where  aluminum  nitride  (AIN)  solid  is  present, 


1  Naval  Weapon*  Center.  Peak  Overpressures  for  Internal  Blast,  by  II.  I  '.  Kinney,  R.  G.  S.  Sewell,  and  K.  J.  Graham. 
China  Lake,  Calif.,  NWC,  June  1979.  (NWC  TP  6089,  publication  UNCLASSIFIED.) 

Naval  Weapons  Center.  Reactive  Metals  in  Internal  Explosions:  The  Combustion  of  Magnesium  in  Air,  by  R.  A. 
Reinhardt.  China  Lake,  Calif..  NWC,  February  1978.  (NWC  I  M  3429,  publication  UNCLASSIFIED.) 

Naval  Weapons  Center.  Adiabatic  Computation  of  Internal  Blast  for  Magnesium-Cased  Charges  in  Air,  by  R.  A. 
Reinhardt.  China  Lake,  Calif.,  NWC,  April  1979.  (NWCTM  3820,  GIDEPE202-1481,  publication  UNCLASSIFIED.) 
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TABLE  1.  Properties  of  the  Fuels. 


Designation 

Name  and  chemical  formula 

^f,298 

(kJ/mole)a 

Relative 

oxygen 

content* 

PETN 

Pentaerythritol  tetranitrate 
C5H8N4°12 

-489.8 

1.33 

NC 

Nitrocellulose,  13.3%  N 
C6H7N2.5°10 

-812.9 

1.05 

HMX 

svm-Cyclotetremethylenetetra- 
nitramine,  C4HgNgOg 

+77.4 

1.00 

RDX 

sym-Cyclotrimethylenetri- 
nitramine,  CjHgNgOg 

+111.47 

1.00 

Pentolite 

50%  PETN,  50%  TNT 
C6.16H6.25N3.41°8.5 

-97.9 

0.92 

Comp  B 

65%  RDX,  35%  TNT 
C1.96H2.53N2.22°2.68 

+30.7 

0.83 

TNT 

2, 4, 6-Trinitrotoluene 

C7H5N3O6 

-11.75 

0.63 

c2h4o 

Ethylene  oxide 

-47.29 

0.25 

"Per  mole  of  formula  indicated.  Data  taken  from  or  computed  from  the  report  of  footnote  1 . 
^Number  of  oxygen  atoms/lnumber  of  carbon  atoms  plus  half  the  number  of  hydrogen  atoms). 


TABLE  2.  Code  of  Condensed  Phases  in  the  Aluminum-Fuel-Air  System. 


Symbol 


Formula,  name,  comments 


A 


AljOj,  aluminum  oxide  (alumina) 

Liquid  above  2315  K,  solid  below  2315  K 


N 

L 

G 


AIN(s),  aluminum  nitride 

Assumed  solid  whenever  produced;  melting 
temperature  is  unknown. 

Alll),  liquid  aluminum  metal 

C(s),  graphite 


C 


AI4C3(s),  aluminum  carbide 


Unstable  above  2700  K.  Below  2700  K,  the  Gibbs 
phase  rule4  requires  that  Alll)  and  C(s)  cannot  both 
be  present.  _ 


cither  as  the  sole  condensed  phase  or  along  with  others.  The  solid  lines  represent  conditions 
where  liquid  or  solid  aluminum  oxide  (AI2O3)  is  the  only  condensed  phase  present. 


4  R  A.  Albert y  and  I-.  Daniels,  Physical  Chemistry,  5th  cd.  (SI  Version).  New  York,  John  Wiley  &  Sons,  1979 
pp.  89-91. 
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BASIS  OF  CALCULATIONS 


From  the  First  Law  of  Thermodynamics  it  follows  that  for  an  adiabatic,  constant- 
volume,  nonflow  process  there  must  be  no  net  change  in  U,  the  internal  energy;  or  the 
internal  energy  of  all  the  products  at  the  adiabatic  flame  temperature,  T,  must  equal  that 
of  the  reactants  at  entering  conditions.  For  each  component  in  the  mixture  of  the  pro¬ 
ducts,  using  the  ideal-gas  assumption  throughout, 


ut  -  u298  +  £  Cyd' 


where  Cy  is  the  heat  capacity  at  constant  volume.  For  U298  the  customary  thermochemical 
convention  may  be  chosen  of  using  the  internal  energy  of  formation  at  298  K.  The  latter  is 
defined  as  the  internal  energy  change  for  the  synthesis  of  the  substance  considered  from  the 
elements  in  their  stable  states  at  the  given  temperature  (298  K).  Of  the  reactant  materials, 
all  but  the  organic  fuel  are  elements  in  the  stable  state  at  298  K,  and  thus  the  initial  internal 
energy  is  just  equal  to  the  internal  energy  of  formation  of  the  fuel,  AUf,  at  298.  Hence,  the 
adiabatic-isochoric  condition  may  be  written: 


S  Ui,T  * 


the  index,  i,  applying  to  all  the  product  species. 

A  temperature  must  be  found  for  which  this  criterion  is  met;  it  may  then  be  identified 
as  the  adiabatic  flame  temperature.  The  pressure  in  the  system  is  then  calculated  from  the 
ideal  gas  law : 

P  =  RT/V  £  n.  (3) 

in  which  nj  is  the  mole  number  for  the  i1*1  product  component,  considering  only  the  gas¬ 
eous  products,  V  the  total  volume  ( 1  cubic  meter)  and  R  the  universal  gas  constant.  The 
overpressure  is  found  by  subtracting  I  bar. 

Thermochemical  data  are  represented  by  a  five-parameter  expression  (of  somewhat 
different  form  than  was  used  previously3): 

UT  =  B,  +  B2T  +  B3T2  +  B4lnT  +  B5/T  (4) 

Data  were  fitted  to  this  equation  using  a  nonlinear  regression  procedure  originally 
designed  by  Dr.  W.  M.  Tolies- 5  Table  3  gives  the  internal  energy  parameters  for  the  products 

5  See  L.P.  Crawford,  “Electron  Paramagnetic  Resonance  of  the  CHjCOj’  Radical  in  Irradiated  Zinc  Acetate  Di¬ 
hydrate,”  M.S.  Theaij,  Naval  Postgraduate  School,  Monterey,  CA.  June  1968,  p.  36. 
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TABLE  3.  Internal  Energies  of  Combustion  Products. 


Expressed  as  a  function  of  T.  the  absolute  temperature  in  kllokelvins  IT=  T/1000);  values  given 
in  joules/mole:  U(T)  =  B1  +  B2T  +  B3T2  +  84  InT  +  B5/T 


Substance 

B1 

B2 

B3 

B4 

B5 

Al 

310931 

11660 

75 

1364 

664 

AIH 

264015 

30151 

272 

-2961 

1040 

AIOH 

-55866 

60797 

-355 

-22628 

-4116 

AUO 

-149141 

50095 

•12 

-572 

1672 

AI  + 

824364 

11395 

70 

2734 

2271 

Ar 

-3718 

12473 

0 

0 

0 

CO 

-84115 

31136 

-25 

-6218 

24 

co2 

-361013 

56789 

67 

-9201 

1764 

H 

213029 

12473 

0 

0 

0 

OH 

166505 

34786 

-88 

20093 

-7842 

h2 

136535 

32890 

313 

-21058 

-9413 

h2o 

50545 

63224 

-395 

-45949 

-16417 

AIO 

844161 

86475 

-3602 

-90479 

-39617 

aio2 

229458 

51802 

226 

3388 

3818 

NO 

105501 

30619 

10 

•4353 

631 

n2 

38785 

31651 

-61 

-7831 

-657 

O 

214445 

9604 

298 

4428 

2050 

o2 

32105 

29476 

662 

-6265 

-2102 

AIO) 

795 

31748 

0 

0 

0 

AIN(S) 

-297116 

37566 

181 

-407 

859 

ai2o3(i) 

1632491 

136549 

0 

0 

0 

AI^Ools) 

1471366 

146423 

1833 

-37420 

12025 

C(sl 

-40329 

22904 

320 

3401 

5429 

CN- 

102390 

31551 

-69 

-7172 

-392 

CN 

691554 

38828 

740 

-37192 

-24658 

c2h 

758476 

66613 

239 

-43882 

-20692 

c2n 

538335 

54371 

-15 

-1052 

2404 

AI4C3(S) 

-89327 

204334 

136 

-31427 

-2969 

HCN 

284944 

62760 

302 

-26216 

-5947 

HNCO 

36816 

83592 

-474 

-28281 

-4750 

HCO 

90276 

56166 

-353 

-19233 

-2346 

ch2o 

54534 

85563 

-609 

-32362 

-1612 

c2h2 

474890 

92019 

-27 

-42593 

10539 

C3 

906056 

48907 

282 

-16185 

-3220 

Substances  are  gaseous  unless  otherwise  specified. 


considered;  data  from  the  JANAF  Tables6  were  used  as  the  primary  source.  For  the  alumi¬ 
num  species  below  2767  K  (the  normal  boiling  point  of  the  metal)  the  thermochemical  data 
were  first  recomputed  to  accommodate  the  universal  choice  of  aluminum  vapor  as  the  ref¬ 
erence  state.  (See  Appendix  A  for  further  discussion  of  this  point.) 


6  National  Bureau  of  Standards,  JANAF  Thermochemical  Tables,  2nd  ed.,  by  D.  R.  Stull  and  H.  Prophet.  Wash¬ 
ington.  D.C.,  NBS,  June  1971. 


6 


NWC  TP  6287 


EQUILIBRIUM  CONSIDERATIONS 


SPECIES  CONSIDERED 

The  equilibrium-stoichiometric  calculations  involved  the  consideration  of  28  gaseous 
species  and  up  to  five  condensed  phases.  The  gaseous  species  are,  rather  arbitrarily,  divided 
into  “major”  and  “minor”  categories: 

Major  (gaseous):  Al,  Al-,0,  Ar,  CO.  CO-,,  H^,  H->0,  Nv  O2 

Minor  (gaseous):  A1H,  AlOH,  Ai+,  H.  OH,  AlO,  AlO,,  NO,  O,  CN*.  CN.  C-,H, 
C,N,  HCN,  HNCO,  CH^O,  C2H^,  C3 

Condensed:  Al(l),  A^O^  (1  or  s),  AlN(s),  Al^C^fs),  C  (graphite). 

Using  the  computational  scheme  described  below,  in  an  abbreviated  version,  the 
following  species  were  eliminated  as  never  reaching  an  appreciable  amount  under  any  of  the 
conditions  to  be  anticipated:  A1HO,  AIHO2,  Ab02,  C(g),  C2.  C2O,  C3O2. 

Ion  formation  was  also  looked  for  in  the  same  way.  It  was  realized  at  once  that  simple 
equilibrium  ionization  to  produce  appreciable  electron  concentrations  would  not  occur, 
even  at  4500  K.  (the  highest  computed  temperature).  But  ions  might  be  formed  in  pairs  by 
the  exchange  of  an  electron  between  neutral  species.  It  was  concluded,  in  this  preliminary 
screening,  that  there  was  the  possibility  of  producing  Al+  and  CN'  at  high  temperatures 
and  under  ideal  stoichiometries. 

A  comment  should  be  made  regarding  the  vaporization  of  the  condensed  phases.  At 
the  temperatures  considered,  liquid  aluminum  is  fairly  volatile,  and  aluminum  vapor  is  thus 
considered  as  one  of  the  important  species.  The  vapor  pressures  of  AlN(s)  to  AlN(g)  or 
of  C(s)  to  either  C(g)  or  C2(g)  are  too  small  to  be  of  significance;  C3(g)  was  included  in  the 
scheme,  though  no  system  was  found  in  which  an  appreciable  amount  occurred.  AI4C3  is 
unstable,  relative  to  the  elements,  above  2700  K.  As  far  as  is  known.  AI2O3  does  not  exist 
as  such  in  the  vapor  phase;  rather,  AI2O3  undergoes  equilibrium  dissociation  at  high  temp¬ 
eratures  into  a  complex  mixture  including  AI2O,  AlO,  O,  etc. 


EQUILIBRIUM  DATA 

For  each  species  there  exists  at  every  temperature  the  equilibrium  constant  of  form¬ 
ation;  i.e. ,  the  equilibrium  constant  for  the  synthesis  of  1  mole  from  the  elements  in  their 
standard  states.  Standard  states  chosen  for  this  purpose  are:  AI(g),  C(s),  H2(g),  N?(g), 
02(g).  Data  were  taken  from  the  JANAF  Tables6  and  fitted  to  the  equation 

log]0K  =  A,  +  A2/(A3  +  T)  +  A4T  (5) 
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using  the  nonlinear  regression  scheme. s  For  aluminum  compounds  below  2700  K,  which  is 
the  boiling  point  of  the  metal,  the  JANAF  data  must  first  be  recalculated  to  account  for  the 
change  in  standard  state  from  liquid  metal  used  in  the  tables  to  the  vapor  as  used  by  us. 
Table  4  gives  the  parameters  for  the  equilibrium-constant  equation  (Equation  5). 

In  application,  the  formation  constants  must  also  be  converted  from  K  (that  is,  in 
terms  of  pressures)  to  Kn  (in  terms  of  mole  numbers).  At  the  beginning  of  any  computation 
at  a  new  temperature,  this  conversion  is  carried  out  from  the  relation 

A  n 

K  =  K  (V/RT)  (6) 

P  n 

where  An  is  the  change  in  number  of  moles  of  gas  for  the  synthesis  reaction.  Then  for  any 
species  the  mole  number  is  given  as  a  known  function  of  the  Kn  of  formation  of  the  mole 
numbers  of  Al,  fH,  N2  and  O2,  and  of  the  activity  of  carbon  (standard  state,  graphite). 
The  activity  concept  is  discussed  in  Appendix  A. 


TABLE  4.  Equilibrium  Constants  of  Formation 
of  Combustion  Products. 


Logarithm  (base  10)  of  the  formation  constant  of  the  indicated  substance,  expressed  as  a  function  of  T,  the 
absolute  temperature  in  kilokelvins  (T  =  T/1000):  LOG(IO)  K  =  Al  +  A2/(A3  +  T)  +  A4r. 


Substance 

Al 

- A2 

A3 

A4 

AIH 

•2.046 

3.096 

-0.061 

0076 

AIOH 

-6.828 

26.915 

0.002 

0.010 

AUO 

-8.582 

40.171 

-0.009 

0.175 

Af+ 

2.255 

-32.796 

0.027 

0.098 

CO 

4.593 

6.108 

0.030 

0.067 

co2 

0.087 

20.642 

0.001 

0024 

H 

3.132 

-12.016 

0019 

0.015 

OH 

0.769 

-1.969 

-0  016 

0.015 

h2o 

-3.056 

13.305 

0014 

0.005 

AIO 

-2.943 

12.989 

:  005 

0.068 

aio2 

•6.148 

26.186 

0  007 

0.010 

NO 

0.710 

-4F0 

0.008 

0.009 

o 

3.497 

-13.4^9 

0.010 

0.006 

akii 

-5.953 

15.658 

0.021 

0.090 

AIN(S) 

-10.936 

31.054 

0.037 

0.117 

ai2o3(i) 

-40.607 

154.393 

0.282 

1.928 

AI203(s) 

-29.436 

118.807 

0.006 

0.506 

CN- 

0.826 

0.005 

-2.418 

0.047 

CN 

5.092 

-22.218 

0.012 

0016 

c2h 

6.529 

-24.269 

0.012 

0.043 

c2n 

7.060 

-29.343 

0.000 

0.066 

Af4C3(Si 

-24.903 

71.143 

0.049 

0.519 

HCN 

1.607 

-6.807 

0.012 

0.008 

HNCO 

-25.740 

96.491 

2.563 

2.862 

HCO 

2.221 

1.209 

0.219 

0.060 

ch2o 

•2.036 

6.770 

0.018 

0.016 

C2H2 

2.655 

■11.332 

0016 

0.000 

_ Si _ 

10.850 

-41.395 

0.016 

0.186 

Substances  are  gaseous  unless  otherwise  specified. 
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EQUILIBRIUM  CALCULATIONS 

For  argon  the  mole  number  is  always  0.4036,  the  number  of  moles  in  1  cubic  meter  of 
standard  air  at  298  K,  1  bar.  For  each  of  the  remaining  five  elements  a  material  balance 
equation  may  be  written,  representing  the  conservation  of  the  number  of  atoms  of  the 
element.  Each  mole  number  is  a  function  of  those  of  the  elements  in  standard  states  and 
of  the  activity  of  carbon;  hence  five  simultaneous  nonlinear  equations  in  five  unknowns 
are  obtained. 

The  number  of  simultaneous  equations  that  must  be  considered  is  reduced  whenever  a 
condensed  phase  is  present,  since  each  condensed  phase  introduces  a  new  restriction.  It  has 
been  found  that,  in  all  cases  which  are  terminal  states  for  an  adiabatic  process  of  the  type 
considered,  there  is  at  least  one  condensed  phase  (most  commonly  AI2O3)  and  sometimes 
more,  up  to  as  many  as  four  in  one  case.  In  brief,  then  it  becomes  necessary  to  solve  a 
system  of  from  two  to  four  simultaneous  nonlinear  equations.  For  this  purpose  the  Newton- 
Raphson  method  was  used,  as  outlined  in  Appendix  B. 

Initial  approximations  for  entering  the  computational  scheme  are  obtained  by  first 
assuming  that  liquid  aluminum  plus  condensed  AI2O3  are  present,  from  which  case  fairly 
simple  iterations,  based  on  chemical  stoichiometric  considerations,  lead  to  solutions  to  the 
equations  without  requiring  the  formal  Newton’s  method  calculations.  Tests  are  then  made 
for  the  actual  presence  of  these  two  phases:  If  both  are  present,  the  calculation  stands;  if 
either  is  absent,  the  Newton’s  method  is  initiated,  using  as  initial  approximations  those 
mole  numbers  found  from  this  first  calculation.  Typically,  three  to  six  iterations  are  made 
in  each  use  of  Newton’s  method.  Further  detail  on  the  computational  scheme  is  intended 
for  a  forthcoming  report. 


RESULTS 


For  all  the  systems  studied,  a  large  number  of  properties  at  the  adiabatic  condition 
were  found.  These  include  adiabatic  temperature  and  overpressure  (which  are  reported  for 
each  case)  and  product  yields,  heat  capacity,  and  isentropic  coefficient  (some  representative 
values  of  which  are  reported).  Total  concentration  (metal  plus  fuel)  was  varied  regularly 
over  the  range  0. 10  to  10.0  kg/m^.  For  the  metal-fuel-air  systems  C/M  was  varied  system¬ 
atically  also  from  0.10  to  10.0. 

Table  2  gave  the  key  to  the  condensed  phases  present,  as  used  throughout  Tables  5 
through  7.  In  Table  5  are  given  overpressures  (bars),  adiabatic  temperatures  (K),  and  pro¬ 
duct  yields  in  mole  %  for  “pure”  aluminum  in  air.*  In  Table  6  adiabatic  temperature  and 
overpressure  are  given  for  each  of  the  eight  chosen  fuels  in  air,  but  in  the  absence  of 
metal.** 


*  Because  of  limitations  in  the  structure  of  the  computer  program,  these  calculations  had  to  be  run  for  A1  + 
0.01%  C.  This  trace  of  carbon  is  quite  trivial  to  the  computed  results. 

**  I  or  this  purpose  a  simplified  program  was  used,  constructed  with  the  aid  of  Kenneth  J.  Graham,  Naval  Weapons 
Center,  in  which  all  reference  to  aluminum  was  omitted. 
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TABLE  5.  Combustion  of  Pure  Aluminum  in  Air: 


Adiabatic  Temperature,  Overpressure,  and  Product  Yields. 


Concen¬ 

tration, 

kg/m  ^ 

Tempera¬ 
ture,  K 

Over¬ 

pressure, 

bars 

Partial  pressure,  bars0 

Product  yields, 
number  of  moles 

At 

ai2o 

AlO 

aio2 

NO 

n2 

(3 

O 

Altl) 

ai2o3(i) 

AIN(s) 

0.1 

2623 

7.2 

0 

0 

0 

0 

0.20 

6.76 

1.13 

0 

— 

195 

— 

0.2 

3813 

10.5 

0.01 

0.01 

0.04 

0.03 

0.50 

9.73 

0.37 

0.63 

- 

3.56 

- 

0.4 

4034 

11.6 

0.53 

0.66 

0.30 

0.03 

0.16 

1047 

0.26 

0.03 

- 

4.16 

- 

1.0 

2987 

9.9 

0.54 

3.24 

0 

0 

0 

7.05 

0 

0 

- 

1.29 

6.17 

2.0 

3158 

10.8 

2.86 

4.45 

0 

0 

0 

4.41 

0 

0 

- 

- 

29.4 

4.0 

3104 

8.4 

3.98 

4.37 

0 

0 

0 

096 

0 

0 

4.34 

— 

55.6 

10.0 

2223 

-0.4 

0.04 

0.48 

0 

0 

0 

0 

0 

0 

293 

49* 

63.0 

a To  the  nearest  0.01  bar.  Balance  is  argon. 
*  Solid  Al203. 


TABLE  6.  Combustion  of  Pure  Fuels  in  Air:  Overpressures 
(With  Adiabatic  Temperatures  in  Parentheses). 


Concentration, 
kg/m 3 

C2H4O 

TNT 

Comp  B 

Pentolite 

RDX 

HMX 

NC 

PETN 

0.1 

7.8 

48 

39 

4.0 

3.3 

3.2 

3.1 

2.8 

(2408) 

(1625) 

(1335) 

(1384) 

(1170) 

(1149) 

(1137) 

(1053) 

0.2 

10.8 

8.2 

6.8 

7.0 

5.9 

5.8 

5.6 

5.1 

(2813) 

(2437) 

(2022) 

(2095) 

(1762) 

(1729) 

(1706) 

(1574) 

0.4 

119 

12  2 

11.3 

11.3 

10.3 

10.1 

9.5 

9.0 

(2275) 

(2972) 

(27551 

(2805) 

(2501) 

(2460) 

(2419) 

(2258) 

1.0 

13.7 

19  5 

20.0 

19.8 

19.3 

19.0 

17.0 

17.2 

(1513) 

(G) 

(3038) 

(3214) 

(3236) 

(3129) 

(3095) 

(2935) 

(2944) 

2  0 

20  2 

29.3 

33.1 

32.7 

32.7 

32.2 

27.2 

28.9 

(1399) 

(G) 

(2821) 

(33631 

(3396) 

(3376) 

(3336) 

(3018) 

(3209) 

4.0 

34.3 

48.0 

58.9 

58.2 

59.3 

58.1 

46.3 

51.8 

(1369) 

(G) 

126171 

(3437) 

(3493) 

(3546) 

(3496) 

(2991) 

(3392) 

10  0 

78.3 

103.1 

136.4 

135.0 

140.0 

136.8 

102.6 

121.5 

(1393) 

(Gla 

(25591 

(G) 

(3499) 

(3582) 

(3715) 

(3651) 

(2929) 

(3578) 

tf(GI  mdicates  solid  carbon  is  present. 


TabU*  7  gives  overpressure  in  bars,  with  the  temperature  in  parentheses  for  the  aluminum- 
luel-air  mixtures.  Lacli  portion  of  the  table  is  devoted  to  a  particular  fuel  and  the  fuels  are 
arranged  from  the  most  oxygen-deficient,  ethylene  oxide,  to  the  most  oxygen-rich,  PETN. 
In  all  the  following,  C  refers  to  the  concentration  of  fuel  (“charge”)  and  M  to  that 
of  metal,  in  kilograms  per  cubic  meter.  C/M  thus  represents  the  charge-to-metal  ratio  and 
C  +  M  the  total  concentration,  in  kilograms  per  cubic  meter.  Again  for  each  table  entry 
is  a  code  indicating  the  condensed  phase  or  phases  present  as  indicated  in  Table  2. 
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TABLE  7.  Overpressure,  in  bars.  Adiabatic  Temperature,  in  Kelvins 
(in  parentheses)  and  Condensed  Phases  for  the 
Combustion  of  Aluminized  Fuels.2 


Fuel:  Ethylene  oxide 


VTM 

0.1 

0.25 

0.5 

i 

2 

4 

10 

0.1 

7.2 

7.3 

7.4 

7.5 

7.6 

7.7 

7.7 

(2592) 

(2563) 

(2534) 

(2499) 

(2467) 

(2443) 

(2425) 

A 

A 

A 

A 

A 

A 

A 

02 

10.5 

10.5 

10.6 

10.7 

103 

103 

103 

(3661) 

(3516) 

(3378) 

(3241) 

(3116) 

(3011) 

(2913) 

A 

A 

A 

A 

A 

A 

A 

0.4 

12.2 

123 

13.6 

14.3 

14.3 

13.7 

123 

(3936) 

(3825) 

(3700) 

(3518) 

(3216) 

(2885) 

(2574) 

A 

A 

A 

A 

A 

A 

A 

1.0 

11.3 

123 

14.7 

17.0 

18.5 

163 

16.3 

(2968) 

(2943) 

(2907) 

(2836) 

(2656) 

(2155) 

(1934) 

AN 

AN 

AN 

AN 

AN 

NG 

ANG 

13.6 

17.0 

20.8 

24.6 

26.7 

25.5 

263 

(3137) 

(3117) 

(3077) 

(2955) 

(2611) 

(2118) 

(1991) 

§H 

N 

N 

N  . 

N 

NG 

NG 

ANG 

1 

13.9 

20.2 

28.6 

37.9 

44.6 

43.4 

45.6 

(3052) 

(2982) 

(2949) 

(2891) 

(2646) 

(2133) 

12045) 

NL 

NLG 

NLG 

NG 

NG 

ANG 

AG 

10.0 

12.2 

29.4 

50.6 

72.2 

94.6 

103.3 

1033 

(2460! 

(2570) 

(2640) 

(2663) 

(2670) 

(2367) 

(2032) 

ALNC 

ALNC 

ALNGC  | 
_ 1 

ALNGC 

ANG 

ANG 

AG 

Fuel:  TNT 


0.1 

7.0 

6.8 

6.6 

6.4 

53 

5.5 

5.1 

(2536) 

(2432) 

(2315) 

(2210) 

(2020) 

(1869) 

(1739) 

A 

A 

A 

A 

A 

A 

A 

03 

10.3 

10.1 

93 

9.4 

9.1 

83 

8.5 

(3666) 

(3479) 

(3276) 

(3062) 

(2865) 

(2703) 

(2565) 

A 

A 

A 

A 

A 

A 

A 

0.4 

12.0 

12.5 

13.0 

13.2 

123 

12.6 

12.3 

(4029) 

(4009) 

(3962) 

(3791) 

(3480) 

(3258) 

(3092) 

A 

A 

A 

A 

A 

A 

A 

1.0 

10.9 

12.0 

13.4 

16.4 

19.4 

203 

20.3 

(2980) 

(2972) 

(2960) 

(3283) 

(3567) 

(3589) 

(3320) 

AN 

AN 

AN 

A 

A 

A 

A 

°See  Table  2  for  condensed  phase*  code. 
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TABLE  7.  Overpressure,  in  ban.  Adiabatic  Temperature,  in  Kelvins 
(in  parentheses)  and  Condensed  Phases  for  the 
Combustion  of  Aluminized  Fuels2  (Contd.). 


Fuel:  TNT 


C+M 

C/M 

0.1 

0.25 

0.5 

1 

2 

4 

10 

2.0 

12.7 

15.0 

17.6 

20.9 

24.2 

26.7 

30.7 

13136) 

(3102) 

(3043) 

(2980) 

(2948) 

(2906) 

(3126) 

N 

N 

N 

AN 

AN 

AN 

A 

4.0 

13.5 

19.7 

26.4 

33.6 

40.4 

45.1 

47.4 

(3154) 

(3196) 

(3220) 

(3190) 

(3118) 

(3015) 

(2857) 

NL 

NL 

N 

N 

N 

N 

N 

10.0 

9.1 

22.3 

40.7 

65.3 

85.5 

98.3 

104.4 

(2529) 

(2829) 

(3090) 

(3235) 

(3234) 

(3137) 

(2909) 

ALNC 

NLG 

NLG 

NLG 

N 

_ 1 _ 1 

NG 

Fuel:  Composition  B 


0.1 

7.0 

6.7 

6.6 

6.0 

5.4 

48 

4.3 

(2515) 

(2385) 

(2315) 

(2086) 

11844) 

(1650) 

(1482) 

A 

A 

A 

A 

A 

A 

A 

10.3 

10.0 

9.6 

9.1 

8.6 

8.0 

79 

(3651) 

(3444) 

(3213) 

(2952) 

(2686) 

(2442) 

(2243) 

A 

A 

A 

A 

A 

A 

A 

0.4 

12.0 

12.5 

13.0 

13.0 

12.5 

12.1 

11.6 

(4042) 

(4032) 

(3976) 

(3738) 

(3390) 

(3139) 

(2930) 

A 

A 

A 

A 

A 

A 

A 

1.0 

11.0 

12,3 

15.5 

18.7 

20.8 

209 

20.4 

(2985) 

(2982) 

(3494) 

(3867) 

(3925) 

(3695) 

(3424) 

AN 

AN 

A 

A 

A 

A 

A 

12.8 

15.1 

18.4 

22.6 

30.4 

339 

339 

(3114) 

(3028) 

(3029) 

(3081) 

(3768) 

(3878) 

(3629) 

N 

AN 

AN 

A 

A 

A 

A 

14.8 

21.6 

27.8 

35.4 

46.2 

578 

599 

(3183) 

(3227) 

(7173) 

(3093) 

(3345) 

(3842) 

(3750) 

NL 

N 

N 

AN 

A 

A 

A 

10.0 

11.1 

29.3 

52.6 

76.2 

96.9 

1252 

1402 

(2613) 

(3104) 

(3285) 

(3255) 

(3212) 

(3715) 

(3858) 

ALN 

NL 

NL 

N 

AN 

A 

A 

Fuel: 

Pentolite 

0.1 

7.0 

6.7 

6.6 

6.1 

Kg 

49 

4.4 

(2518) 

(2393) 

(2315) 

(2106) 

(1686) 

(1524) 

A 

_ * _ 1 

A 

A 

mm 

A 

A 

2See  Table  2  for  condensed  phases  code. 
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TABLE  7.  Overpressure,  in  bars,  Adiabatic  Temperature,  in  Kelvins 
(in  parentheses)  and  Condensed  Phases  for  the 
Combustion  of  Aluminized  Fuels?  (Contd.). 


Fuel:  Pentolite 


C+M 

C/M 

0.1 

0.25 

0.5 

1 

2 

4 

10 

0.2 

10.3 

9.9 

9.6 

9.1 

8.6 

8.1 

7.6 

<3651) 

(3445) 

(3218) 

(2968) 

(2716) 

(2489) 

(2304) 

A 

A 

A 

A 

A 

A 

A 

0.4 

12.0 

12.5 

12.9 

12.9 

12.4 

12.0 

11.6 

(4045) 

(4037) 

(3979) 

(3735) 

(3394) 

(3156) 

(2962) 

A 

A 

A 

A 

A 

A 

A 

1.0 

11.0 

12.3 

16.7 

18.7 

20.6 

20.6 

202 

(2984) 

(2981) 

(3570) 

(3918) 

(3946) 

(3698) 

(3434) 

AN 

AN 

A 

A 

A 

A 

A 

2.0 

126 

15.1 

18.3 

23.5 

30.6 

33.5 

33.4 

(3108) 

(3025) 

(3025) 

(3268) 

(3873) 

(3918) 

(3649) 

N 

AN 

AN 

A 

A 

A 

A 

14.6 

21.4 

27.4 

35.1 

476 

57.6 

586 

(3166) 

(3205) 

(3138) 

(3084) 

(3573) 

(3957) 

(3793) 

NL 

N 

N 

AN 

A 

A 

A 

10.0 

10.4 

27.7 

50.7 

74.3 

95.8 

127  6 

138.4 

(2568) 

(2953) 

(3193) 

(3178) 

(3197) 

(3929) 

(3934) 

ALN 

NL 

NL 

AN 

AN 

A 

A 

Fuel:  RDX 


0.1 

6.9 

6.6 

6.5 

5.8 

5.1 

4.4 

36 

(2504) 

(2360) 

(2272) 

(2016) 

(1745) 

(1527) 

(1336) 

A 

A 

A 

A 

A 

A 

A 

0.2 

10.2 

9.9 

9.5 

8.9 

8.2 

7.6 

6.7 

(3643) 

(3424) 

(3175) 

(2881) 

(2564) 

(2311) 

(2024) 

A 

A 

A 

A 

A 

A 

A 

0.4 

12.0 

12.5 

13.0 

12.9 

12.2 

11.7 

106 

(40491 

(4042) 

(3976) 

(3700) 

(3322) 

(3031) 

(2758) 

A 

A 

A 

A 

A 

A 

A 

1.0 

11.1 

12.9 

16.4 

19.2 

20.8 

20.5 

196 

(2987) 

(3141) 

(3734) 

(4008) 

(3962) 

(3661) 

(3369) 

AN 

A 

A 

A 

A 

A 

A 

126 

15.4 

186 

26.2 

32.0 

34.0 

33.5 

(3100) 

(3037) 

(3043) 

(3672) 

(4051) 

(3955) 

(3646) 

N 

AN 

AN 

A  1 

A 

A 

A 

JSee  Table  2  for  condensed  phases  code. 
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TABLE  7.  Overpressure,  in  bars,  Adiabatic  Temperature,  in  Kelvins 
(in  parentheses)  and  Condensed  Phases  for  the 
Combustion  of  Aluminized  Fuels®  (Contd.). 


Fuel:  RDX 


C/M 

0.1 

0.25 

0.5 

1 

2 

4 

10 

4.0 

15.5 

22.1 

28.2 

37.0 

53.1 

60.0 

59.9 

(3197) 

(3217) 

(3108) 

(3123) 

(3993) 

(4119) 

(3843) 

NL 

N 

AN 

AN 

A 

A 

A 

10.0 

12.4 

33.2 

57.8 

79.6 

115.2 

138.0 

142.6 

(2745) 

(3189) 

(3328) 

(3253) 

(3904) 

(4272) 

(4052) 

_ i 

NL 

NL 

N 

AN 

A 

A 

A 

Fuel:  HMX 


0.1 

65 

6.6 

6.5 

59 

5.0 

4.3 

3.7 

(2356) 

(2266) 

(2007) 

(1732) 

(1511) 

(1318) 

A 

A 

A 

A 

A 

A 

A 

0.2 

10.2 

95 

9.5 

8.9 

8.1 

7.5 

6.6 

(3641) 

(3420) 

(3168) 

(2870) 

(2547) 

(2288) 

(1995) 

A 

A 

A 

A 

A 

A 

A 

0.4 

12.0 

12.5 

125 

128 

12.2 

11.6 

109 

(4049) 

(4041) 

(3973) 

13689) 

(3306) 

(3010) 

(2728) 

A 

A 

A 

A 

A 

A 

A 

1.0 

11.1 

128 

16.3 

19.1 

20.6 

20.3 

19.7 

(2987) 

(3124) 

(3720) 

(3999) 

13948) 

(3637) 

(3340) 

AN 

A 

A 

A 

A 

A 

A 

2.0 

12.9 

15.3 

18.7 

255 

31.7 

33.6 

33.0 

(3099) 

(3036) 

(3041) 

(3641) 

(4033) 

(3929) 

(3610) 

N 

AN 

AN 

A 

A 

A 

A 

4.0 

15.4 

215 

28.0 

36.7 

52.4 

595 

589 

(3196) 

(3214) 

(3106) 

(3120) 

(3961) 

(4091) 

(3800) 

NL 

N 

AN 

AN 

A 

A 

A 

10.0 

12.3 

265 

46.5 

78.8 

113.1 

136.0 

1398 

(2733) 

(2985) 

(3151) 

(3249) 

(3850) 

(4235) 

(3997) 

NL 

NL 

NL 

_ 1 

AN 

A 

A 

A 

Fuel:  Nitrocellulose 

0.1 

65 

6.6 

6.4 

1 

5.7 

4.3 

3.6 

(2501) 

(2354) 

(2261) 

(200 1) 

(1501) 

(1308) 

A 

A 

A 

A 

A 

A 

0.2 

10.2 

9.8 

9.3 

8.7 

7.3 

6.4 

(3630) 

(3400) 

(3144) 

(2848) 

(2270) 

(1974) 

A 

A 

A 

A 

■of 

A 

A 

®See  Table  2  for  condensed  phases  code. 
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TABLE  7.  Overpressure,  in  bers,  Adiebetic  Temperature,  in  Kelvins 
(in  parentheses)  and  Condensed  Phases  for  the 
Combustion  of  Aluminized  Fuels?  (Contd.). 


Fuel:  Nitrocellulose 


C+M 

C/M 

0.1 

0.26 

0.5 

1 

1 

2 

4 

10 

0.4 

12.0 

12.4 

12.7 

12.3 

1  1 

11.6 

11.0 

109 

(4044) 

(4030) 

(3944) 

(3609) 

(3226) 

(2944) 

(2679) 

A 

A 

A 

A 

A 

A 

A 

1.0 

10.9 

12.2 

15.5 

18.1 

19.1 

18.3 

17.6 

(2983) 

(3010) 

(3623) 

(3932) 

(3816) 

(3449) 

(3163) 

AN 

A 

A 

A 

A 

A 

A 

2.0 

12.5 

14.7 

17.7 

23.3 

28.9 

295 

28.3 

(3077) 

(3020) 

(3018)- — 

(3416) 

(3897) 

(3672) 

(3309) 

N 

AN 

AN 

A 

A 

A 

A 

14.5 

20.4 

25.8 

33.4 

46.2 

51 2 

48.7 

(3148) 

(3150) 

(3062) 

(3068) 

(3718) 

(3801) 

(3375) 

NL 

N 

AN 

AN 

A 

A 

A 

10.0 

9.8 

26.1 

48.1 

69.8 

94.6 

1159 

111.4 

(2569) 

(2760) 

(3021) 

(3153) 

(3436) 

(3899) 

(3424) 

ALN 

ALN 

NL 

AN 

A 

A 

A 

0.1 

6.9 

6.6 

6.4 

5.6 

4.8 

4.1 

3.4 

(2495) 

'2342) 

(2241) 

(1967) 

(1675) 

(1440) 

(1234) 

A 

A 

A 

A 

A 

A 

A 

0.2 

10.2 

99 

9.3 

8.6 

79 

7.0 

6.1 

(3630) 

(3397) 

(3133) 

(2817) 

(2467) 

(2182) 

(1863) 

A 

A 

A 

A 

A 

A 

A 

0.4 

12.0 

12.4 

12.7 

12.3 

11.6 

109 

10.0 

(4058) 

(4052) 

(3962) 

(3618) 

(3219) 

(2902) 

(2580) 

A 

A 

A 

A 

A 

A 

A 

1.0 

11.1 

13.6 

16.6 

185 

19.5 

189 

18.0 

(2987) 

(3401) 

(3920) 

(4098) 

(3895) 

(3521) 

(3215) 

AN 

A 

A 

A 

A 

A 

A 

2.0 

12.7 

15.3 

189 

265 

30.6 

305 

299 

(3066) 

(3036) 

(3089) 

(4028) 

(4162) 

(3844) 

(3493) 

N 

AN 

A 

A 

A 

A 

A 

4.0 

15.5 

21.4 

279 

41.4 

519 

54.4 

525 

(3172) 

(3144) 

(3103) 

(3803) 

(4290) 

(4090) 

(3701) 

NL 

N 

AN 

A 

A 

A 

A 

flSee  Table  2  for  condensed  phases  code. 


-  J 


t 


Fuel:  PETN 


C+M 

r 

0.1 

0.25 

0.6 

1 

2 

4 

10 

|  10.0 

11.8 

32.4 

55.1 

81.8 

116.4 

126.3 

124.6 

(2656) 

(3016) 

(3171) 

(3464) 

(4443) 

(4353) 

(3929) 

ALN 

NL 

AN 

A 

A 

A 

A 

aSee  Table  2  for  condensed  chases  code. 


In  Table  8  is  given  the  set  of  product  yields  for  the  TNT-aluminum-air  system  at 
C/M  =  I  0,  chosen  as  representative.  Mole  %  in  the  vapor  is  given  for  all  components  that 
reach  at  least  0.1%  at  some  point.  Numbers  of  moles  of  condensed  phases  are  also  provided 
Also  given  are  Cy,  the  specific  heat  capacity  at  constant  volume,  in  kJ/kg/K  and  an  isen- 
tropic  parameter  A,  defined  such  that  TV^‘1  is  constant  for  an  isentropic  process  of  small 
displacement.  Computation  and  significance  of  Cy  and  A  are  discussed  in  Appendix  C.  It 
should  be  added  that  product-yield  data  and  values  of  Cy  and  A  are  available  for  all  the  sys¬ 
tems  run  but  have  not  been  included  in  this  report  due  to  the  bulk  of  data  involved. 


DISCUSSION 


PURE  ALUMINUM  IN  AIR 

Table  5,  as  remarked,  gives  the  data  for  aluminum-air  mixtures.  It  is  noted  that 
temperature  rs  concentration  shows  a  maximum  at  0.4  kg/m^.  This  corresponds  reason¬ 
ably  well  with  the  notion  that  at  0.305  kg/m^,  aluminum  and  oxygen  are  present  in  stoi- 
chiometrically  equivalent  amounts  to  form  AI2O3.  Beyond  this  point,  to  an  increasing 
extent,  Al203(l)  is  being  converted  to  A^Ofg)  by  the  endothermic  (energy-absorbing) 
reaction 


A1,03(1)  +  4  Al  —  3  Al,0(g)  (7) 

The  drop  in  temperature  is  arrested  when  aluminum  nitride  begins  to  form.  This  exothermic 
process. 


Ai  +  1/2  Ny— -  AlN(s) 


(8) 
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TABLE  8.  Combustion  of  TNT-Aluminum  in  Air  at  C/M  -  1.0:  Temperature, 
Prenure,  Product  Yield  and  Itentropic  Parameter*. 


C  +  M,  kg/m3 


0.1 

0.2 

0.4 

1.0 

2.0 

4.0 

10 

T,  K 

2210 

3062 

3791 

3283 

2980 

3190 

3235 

Overpressure, 

bars 

6.4 

9.5 

13.2 

16.4 

30.1 

33.6 

65.3 

Al 

_ 

— 

0.1 

2.4 

2.2 

8.9 

9.4 

AIH 

- 

- 

- 

0.1 

0.2 

1.0 

1.5 

AIOH 

- 

- 

- 

0.1 

0.1 

0 

0 

ai2o 

- 

- 

0.1 

8.4 

12.7 

8.1 

8.4 

CO 

2.4 

12,2 

24.1 

349 

45.4 

52.0 

co2 

3.8 

5.1 

1.5 

0.1 

0 

0 

0 

H 

- 

0.2 

3.1 

2.3 

1.1 

1.8 

1.6 

£ 

OH 

0.1 

1.1 

1.7 

0 

0 

0 

0 

C 

h2 

- 

0.1 

1.2 

7.2 

11.7 

14.6 

17.0 

i 

h2o 

1.3 

1.9 

1.2 

0.1 

- 

- 

- 

| 

AIO 

- 

- 

0.3 

0.1 

- 

- 

- 

aio2 

- 

- 

0.1 

- 

- 

- 

- 

NO 

1.1 

31  j 

2.2 

- 

- 

- 

- 

n2 

78  4 

76.8 

71.5 

54.5 

36.6 

18.2 

4.6 

O 

0 

1.1 

2.7 

- 

- 

- 

°2 

14.1 

7.2 

1.1 

- 

- 

- 

- 

HCN 

- 

- 

- 

0.1 

1.4 

2.6 

C2H2 

- 

-  i 

_ 

- 

- 

0.2 

2.3 

Moles  of 

99.7* 

AINIsi 

- 

- 

- 

11.6 

40.0 

Moles  of 

1 

Al,0,tl) 

0.9^ 

1.9 

3.5 

3.1 

0.4 

- 

- 

cv. 

J/K/fl 

1.18 

2.12 

3.72 

1.69 

25.8 

12.2 

10.3 

1.23 

1.14 

1.10 

1.17 

1.07 

1.07 

1.06 

"Balance  it  argon. 

hPlus  17.2  molet  AKU  and  6.0  molej  C(t). 

"Solid  A^Oj. 

'^The  itentropic  coefficient  defined  in  Appendix  C. 


is  not  possible  until  the  free  oxygen  has  been  consumed  in  the  formation  of  aluminum 
oxide.  The  final  drop  in  temperature  at  the  high-concentration  end  is  the  result  of  the  ap¬ 
pearance  of  liquid  aluminum  metal  in  large  quantities,  which  serves  as  an  energy  sink  due 
solely  to  its  heat  capacity. 

The  curve  of  total  pressure  vs  concentration,  as  shown  subsequently,  is  labeled 
“C/M  =  0.”  At  low  concentrations  the  pressures  reflect  the  corresponding  temperature 
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effects,  but  beyond  the  equivalance  point  (0.3  kg/m^)  the  pressure  does  not  change  as  much 
as  the  temperatures  might  suggest,  largely  a  result  of  the  increased  number  of  moles  of  gaseous 
AbO.  At  the  high-concentration  extreme  the  surprising  circumstance  is  noted  that  the  total 
pressure  drops  to  less  than  atmospheric  (resulting  in  a  negative  overpressure)  At  this  point 
only  about  1%  of  the  material  present  is  in  the  gas  phase. 


TEMPERATURE  AND  PRESSURE  IN  ALUMINUM-FUEL-AIR  MIXTURES* 

Figure  1  shows  temperature  vs  total  concentration  (C  +  M)  for  two  representative 
C/M  ratios  for  TNT  (an  oxygen-deficient  explosive)  and  PETN  (an  oxygen-rich  explosive). 
For  the  most  part,  the  maximum  in  temperature  occurs  fairly  close  to  the  equivalence 
point  for  the  formation  of  AI2O3  +  CO  +  H2  +  N2  (marked  E  on  the  curves)  and  for 
approximately  the  same  reason  as  with  aluminum  itself.  But  for  PETN  of  low  aluminum 
content  large  C/M  (and  the  same  is  observed  for  other  oxygen-rich  fuels)-no  maximum 
is  observed  even  though  the  systems  become  oxygen-deficient  at  high  concentrations.  There 
is  sufficient  oxygen  in  the  fuel  so  that  as  C  +  M  increases  the  principal  effect  is  the 
exothermic  production  of  AI2O  from  its  elements.  This  results  in  a  slowing  of  the  temp¬ 
erature  rise,  since  the  process  is  less  exothermic  than  the  formation  of  AI2O3,  but  sufficient 
so  as  to  overcome  the  energy-wasting  production  of  H2.  The  lower  right-hand  part  of  the 
figure,  shown  by  dotted  lines,  is  that  where  AIN  is  present. 

The  pronounced  dependence  of  total  pressure  on  total  concentration  is  shown  in 
Figure  2,  for  TNT  at  two  different  C/M  ratios.  These  curves  are  explainable  on  the  basis 


FIGURE  1.  Adiabatic  Tamparatura  vt  Total  Concentration  for  TNT  and  PETN 
With  Aluminum  and  Air. 


*  Note  that  in  figures  1  through  5,  dotted  lines  arc  used  to  indicate  the  presence  of  AIN  solid. 
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FIGURE  2.  Total  Franura  t*  Total  Conoantration  of  Pura  TNT,  Pure  Aluminum 
and  TNT-Aluminum  Mixture*  in  Air. 

of  the  temperature  effects  just  discussed,  taking  into  account  the  increased  production 
of  gaseous  products  (CO,  H2,  and  AI2O)  at  higher  concentrations.  Table  8  shows  this 
effect,  especially  if  one  multiplies  the  mole  %s  by  the  total  mole  numbers.  It  is  interesting 
to  note,  by  comparison  with  the  temperature  effects  in  Figure  1.  that  at  high  concentra¬ 
tion  the  overpressure  continues  to  rise  rapidly  at  C/M  =  2  but  not  for  C/M  =  0.25,  even 
though  the  temperatures  behave  similarly  in  this  region  (where  AIN  is  the  condensed  phase 
present).  For  comparison  purposes  the  curves  for  pure  TNT  (C/M  =  °°)  and  also  for  pure 
A1  (C/M  =  0)  which  was  discussed  earlier  are  included  in  the  figure. 

The  effects  of  C/M  ratio  on  temperature  and  overpressure  are  shown  in  Figure  3  for 
PFTN  at  two  different  total  concentrations.  Curves  for  other  fuels  will  differ  only  in  details. 
The  decrease  in  temperature  as  C/M  increases  at  low  concentration  (C  +  M  =  0.2  kg/m^)  is 
characteristic  of  the  oxygen-rich  regime.  Aluminum  has  a  higher  heating  value  ( ca  31  kJ/g) 
than  any  of  the  fuels  considered  here  ;  thus,  replacing  aluminum  by  fuel,  as  C/M  is  increased, 
would  result  in  a  lowering  of  temperature.  (A  hydrocarbon,  however,  would  have  shown  the 
reverse  trend.)  At  the  higher  concentration  the  mixture  is  oxygen-deficient  for  small  values 
of  C/M;  increasing  the  content  of  PETN  (or  other  explosive  fuel)  improves  the  oxygen 
balance  and  the  temperature  is  seen  to  rise  to  a  maximum;  this  maximum  in  many  cases,  but 
not  all,  corresponds  closely  to  the  equivalence  point  (AI2O3  +  CO).  The  pressure  effects 
seen  in  Figure  3  parallel  at  least  qualitatively  those  of  temperature,  though  the  pressure 
maximum  at  C  +  M  =  2.0  kg/m 3  is  very  shallow. 

The  effect  of  different  fuels  is  given  for  a  variety  of  total  concentrations,  all  at  C/M  =  1 .0, 
for  the  temperature  in  Figure  4  and  for  overpressure  in  Figure  5.  The  ordinate  for  these 
plots  is  the  relative  oxygen  content  (ROC’)  defined  previously3  as  the  number  of  oxygen 
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FIGURE  3.  Adiabatic  Temperature  and  Total  Pratwra  vt  Charge-  to- Metal  Ratio  for 
PETN  With  Aluminum  and  Air. 

atoms  divided  by  the  sum  of  the  number  of  carbon  atoms  plus  half  the  number  of  hydrogen 
atoms.  ROC  =  1.0  corresponds  to  a  material  just  balanced  for  the  formation  of  CO  +  H^O. 
The  complex  appearance  of  the  temperature  plots  of  Figure  4  is  in  part  explainable.  At  low 
concentration  (C  +  M  =  0.2  kg/m-*)  all  the  systems  are  oxygen-rich.  The  additional  oxygen 
supplied  by  the  fuels  of  higher  ROC  is  thus  a  liability ;  which  is  to  say  that  the  heating  value 
per  unit  mass  decreases  as  the  oxygen  content  increases.  Hven  at  C  +  M  =  0.4  kg/m^  the 
Csl^O-AI-air  system  is  somewhat  oxygen-deficient,  resulting  in  the  drop  in  temperature 
from  TNT  to  ethylene  oxide  along  this  curve.  At  the  higher  concentrations  the  systems  are 
mostly  oxygen-deficient,  and  increasing  the  oxygen  content  alleviates  this  somewhat,  re¬ 
sulting  in  generally  increasing  temperature  with  ROC.  The  irregularity  between  nitrocellulose 
and  RDX  (HMX  is  virtually  identical  in  behavior  to  RDX)  is  presumably  the  result  of  the 
positive  energy  of  formation  of  RDX  (111  kJ/mole)  contrasted  to  the  considerable  negative 
energy  of  formation  of  nitrocellulose  (-813  kJ/mole). 

I'lie  pressure  effects,  in  Figure  5,  are  much  less  striking,  albeit  parallel.  In  particular 
can  be  noted  the  lowered  pressure  for  oxygen-rich  fuels  at  low  concentration  and  the 
enhanced  pressure  for  PFTN,  relative  to  all  the  others,  at  high  concentration. 


CONDENSED  PHASES 


If  is  >|uitc  evident,  from  an  examination  of  the  tables  and  figures,  that  the  nature  of 
the  condensed  phases  present  can  make  profound  alterations  in  the  way  in  which  adiabatic 
temperature  and  overpressure  vary  with  concentration  or  C/M  ratio.  For  this  reason  it  is  of 
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FIGURE  4.  Adiabatic  Temperature  vs  Relative  Oxygen  Content  for  Combustion  of  Eight 
Fuels  in  the  Presence  of  Aluminum  and  Air. 

considerable  interest  to  observe  which  phases  are  present  under  various  conditions  and  to 
attempt  to  consider  how  temperature  and  pressure  might  thereby  be  affected. 

In  Figure  6  is  shown  approximately  the  regions  where  various  condensed  phases  exist 
under  adiabatic  conditions  for  the  TNT-aluminum-airsystem.*  Other  fuels  will  resemble  TNT 
qualitatively,  but  the  regions  will  be  skewed  at  high  C/M  values  according  to  the  oxygen 
content.  For  example,  AIN  is  restricted  to  the  lower  right-hand  20%  of  the  diagram  for  PETN. 

It  may  be  considered  that,  as  C  +  M  is  increased,  the  effect  is  that  of  adding  the  ele¬ 
ments  as  Al(l).  C(s).  Hi,  Ni.  and  Or;  any  oxygen-deficiency  will  be  enhanced  as  C  +  M  is 
increased,  thus  making  it  easier  for  AIN  to  be  formed.  Exothermic  processes  will  tend  to 
raise  the  temperature,  and  thus  the  pressure.  Endothermic  processes  clearly  lower  the  tem¬ 
perature.  but  so  also  do  energetically  neutral  ones,  such  as  the  production  of  Hi  as  a  product 
(as  occurs  in  the  AIN  regions),  because,  as  remarked  earlier,  energy  must  be  supplied  to 
bring  the  Hi  to  the  prevailing  temperature. 

*  A  more  conventional  isothermal  phase-equilibrium  diagram  can  also  be  constructed.  Actual  calculation  has  shown 
that  the  isothermal  diagram  at  3000  K  differs  only  in  detail  from  Figure  6  notably  in  reduction  of  the  rone  where  both 
AIN  and  AlyOj  are  to  be  found. 
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FIGURE  5.  Total  Pressure  vs  Relative  Oxygen  Content  for  Combustion  of  Eight  Fuels 
in  the  Presence  of  Aluminum  and  Air. 

The  effect  on  pressure  is  yet  more  complex,  since  it  depends  not  only  on  temperature 
hut  also  on  the  total  number  of  moles  of  gas.  Examined  below  are  several  important  re¬ 
actions  which  account  for  phase  changes  when  C  +  M  is  increased.  It  may  be  noted  that  for 
each  case,  because  of  conflicting  effects,  it  is  not  possible  to  predict  the  pressure  effect. 

1 .  Simple  production  of  aluminum  nitride: 

Al(l)  +  l/2N2(g)  —  AlN(s) 

The  reaction  is  strongly  exothermic  (though  less  so  than  the  production  of  AI2O3)  and 
thus  increases  the  temperature.  The  number  of  moles  of  gas  decreases. 

2.  Production  of  aluminum  nitride  at  the  expense  of  aluminum  oxide: 

Al-,03(l)  +  N,(g)  +  3C(s)  2AlN(s)  +  3CO(g) 

This  reaction  is  endothermic  but  increases  the  number  of  moles  of  gas. 


22 


NWC  TP  6287 


FIGURE  6.  Condentad  Pham  in  the  Combustion  Products  of  TNT  and  Aluminum  in  Air. 

3.  Disappearance  of  aluminum  oxide,  without  concurrent  formation  of  aluminum 
nitride: 


Al203(!)  +  2C(s)  —  Al20(g)  +  2CO(g) 

A1-,03(I)  +  4AKD— —  3Al?0(g) 

Both  reactions  are  endothermic  but  with  a  marked  increase  in  the  number  of  moles  of  gas. 

A  curious  feature  of  Figure  6  is  the  “reappearance”  of  AbC>3  in  the  extreme  lower 
right-hand  comer.  It  is  instructive  to  examine  the  point  (C/M  =  0. 1 ,  C  +  M  =  10  kg/m^)  for 
all  the  fuels  run,  as  shown  in  Table  9.  It  appears  that  the  presence  or  absence  of  AI2O3O) 
is  governed  entirely  by  the  equilibrium 

A1203(1)  +  4A1(1)  —  3Al2Q(g) 
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TABL£  9.  Properties  of  Aluminum-Fuel-Air  Mixtures 
*t  C/M  =  0.1,  C  +  M  =  10  kfl /m3. 


Fuel 

T 

Condensed  phases 

„  a  “ 

•0  aC  • 
moles 

PETN 

1 

2656 

37 

NC 

2569 

Al(l).  AIN(s),  AljOjtl) 

31 

Pentolite 

2567 

25 

Comp  B 

2613 

23 

RDX 

2745 

AI(l),AIN(s) 

29 

HMX 

2733 

29 

TNT 

2529 

Aid)  ,  AIN(s),  Al203  d),  AI4C3(sl 

13 

c2h4o 

2459 

A 

aSee  text. 


The  formation-constant  data  are  also  used  to  compute  the  number  of  moles  of  AI2O  gas  in 
equilibrium  with  the  two  liquids  at  several  temperatures:  2500  K  (13.4  moles);  2600  K 
(22.2  moles);  2700  K  (35.0  moles).  Thermodynamic  principles  require  that  the  actual  num¬ 
ber  of  moles  of  AhO  may  not  exceed  the  equilibrium  value;  if  smaller  than  the  equilibrium 
value,  then  equilibrium  is  not  established  and.  for  this  case  of  large  excess  of  aluminum,  the 
AI2O3  must  disappear. 

The  oxygen  present  in  these  systems  is.  for  all  but  the  last  two  cases,  converted  first 
to  CO.  then  to  AhO.  First  assuming  no  AI2O3  is  formed,  the  mole  number  of 
AhO  =  aQ  -  a(7  (moles  of  oxygen  atoms  less  moles  of  carbon  atoms  in  the  mixture).  If 
aO  -  ac  <  (AhO)eq  (the  equilibrium  value  at  the  adiabatic  temperature),  then  AI2O3  will 
not  be  observed.  Such  is  the  case  with  HMX  and  RDX  at  ca  2700  K.  PETN  is  so  oxygen-rich 
that  even  at  2700  K.  AI2O3  would  form;  with  nitrocellulose,  pentolite.  and  Comp  B.  the 
adiabatic  temperature  is  significantly  lower,  so  that  again  aQ  -  ac  exceeds  (Al20)eq  The 
cases  of  RDX  and  HMX  are  thus  another  result  of  the  considerably  positive  energy  of 
formation. 

The  very  oxygen-deficient  TNT  and  C2H4O  show  another  feature  altogether.  Here  the 
temperature  is  so  low  and  the  carbon  content  so  high  that  it  is  possible  for  AI4C3  to  be 
produced.  As  a  result  there  is  made  available  considerable  oxygen  that  would  otherwise  have 
been  used  in  formation  of  CO.  and  once  again  AI2O3  can  be  produced.  (In  the  isothermal 
diagram  for  TNT  at  3000  K.  cited  in  the  footnote  on  page  21,  Ah03  is  not  seen.) 


COMMENTS  ON  THE  APPROXIMATIONS 


The  assumption  that  the  process  is  adiabatic  implies  negligible  energy  losses  due  to 
radiation,  conduction,  and  convection  from  the  moment  of  initiation  to  the  establishment 
of  thermodynamic  equilibrium.  Of  these,  the  losses  due  to  radiation  seem  the  more  likely. 
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because  of  the  very  high  temperatures  involved.  A  recent  study  by  LCDR  D.  E.  Smith7 
indicates,  however,  that  where  detonation  has  occurred  the  energy  loss  due  to  radiation  will 
amount  to  less  than  1/2%  of  the  energy  liberated  isothermally.  The  model  assumes  black- 
body  radiation,  an  assumption  which  is  supported  in  the  treatment,  and  considers  the  time 
for  radiation  to  be  no  greater  than  20  times  the  time  required  for  arrival  of  the  detonation 
front  to  the  surface  of  the  container. 

Where  deflagration  is  the  sole  combustion  process,  radiation  effects  may  be  expected 
to  be  considerable.  Some  consideration  of  this  point  has  been  made  in  the  study  of  non- 
metallized  fuels  by  Kinney.  Sewell  and  Graham.1 

The  ideal  gas  assumption  is  very  good  where  the  permanent  gases  are  concerned:  CO 
and  Ht  are  at  temperatures  one  to  two  orders  of  magnitude  greater  that  the  critical  tem¬ 
peratures,  and  maximum  pressures  are  at  most  a  few  times  the  critical  value.  Since  metallic 
aluminum  is  in  the  two-phase  region,  appreciable  deviation  from  ideality  must  be 
considered;  its  critical  temperature  seems  not  to  be  known,  but  a  rough  guess  of  5000  K 
would  indicate  a  severe  departure  from  ideality  in  these  systems. 

It  is  intended  to  address  the  question  of  corrections  for  nonideality  at  a  future  stage 
of  the  project.  The  effects  of  pressure  on  the  apparent  equilibrium  constants  as  well  as  the 
effects  on  computation  of  the  final  pressure  must  be  taken  into  account. 


O.  li.  Smith,  "Attenuation  Effects  of  Thermal  Radiation  on  Internal  Blast  Overpressure,”  M.S.  Thesis,  Naval 
Postgraduate  School,  December  1979. 
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Appendix  A 

STANDARD  STATES,  ACTIVITIES,  AND  FORMATION  CONSTANTS 


The  activity  of  any  chemical  species  in  a  system  may  be  defined6  as 

a.  =  f./f°  S  P-P°  (A- 1) 

i  11  11 

where  the  fs  are  the  fugacities  and  the  approximation  to  partial  pressure  assumes  ideal-gas 
behavior,  the  superscript  °  refers  to  a  standard  state  which  may  be  chosen  arbitrarily. 
Most  commonly  the  standard  state  is  that  of  unit  fugacity  for  gases  (whence  activity  is  ap¬ 
proximated  by  the  partial  pressure)  and  for  the  pure  phase  where  condensed  phases  are 
concerned.  In  addition. reference  states  are  commonly  used  in  thermochemistry  and  thermo¬ 
dynamics,  usually  selecting  the  stable  form  of  the  element. 

On  this  basis,  we  can  write  for  any  chemical  species.  A;  =  AlaC^Hj-N^Oe  the  thermo¬ 
dynamic  equilibrium  constant  of  formation  (from  the  elements  in  their  reference  states) 
as 


Kr  =  a.  /(«  A,aa„ba 


c/2  d/2  e/2 


Arc  h2  n2 


O, 


(A-2) 


bearing  in  mind  that  the  reference  states  for  hydrogen,  nitrogen  and  oxygen  are  the  dia¬ 
tomic  gases. 

In  thermodynamic  tables,  such  as  the  JANAF  Tables,  the  reference  state  is  commonly 
taken  as  the  phase  stable  at  the  temperature  of  concern,  and  the  standard  state  will  change 
accordingly.  Thus,  for  metallic  aluminum  the  standard  state  is  the  pure  liquid  between  the 
melting  point  (933  K)  and  the  normal  boiling  point  (2767  K)  and  is  the  state  of  unit  fu¬ 
gacity  above  that  temperature.  This  change  in  standard  state  leads  to  a  discontinuity  in  the 
slope  of  the  equilibrium  constants  and  would  require  a  different  empirical  equation  to  fit 
the  data  below  2767  K  than  above.  It  was  thus  felt  to  be  far  more  convenient  to  select 
aluminum  vapor  at  unit  fugacity  as  the  standard  state  throughout.  This  choice  did,  however, 
require  recalculation  of  some  of  the  JANAF  formation-constant  data  before  fitting.  Since 
aluminum  vapor  was  frequently  an  important  component  of  the  gaseous  products,  the 
choice  was  quite  a  reasonable  one. 

In  the  case  of  carbon,  however,  it  scarcely  seemed  appropriate  to  select  any  of  the 
gaseous  elemental  species  (C.  C2,  orC3)  as  the  standard  state,  since  this  would  require  re¬ 
computing  all  of  the  data  on  carbon  compounds;  and.  moreover,  since  the  activity  of  car¬ 
bon  becomes  very  small  at  low  temperatures  in  the  presence  of  excess  oxygen,  it  is  entirely 
possible  that  the  amount  of  gaseous  carbon  would  be  in  error  on  account  of  computational 
underflow.  (The  HP  9845  cannot  treat  numbers  less  than  10'99.)  Graphite  was  therefore 
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selected  as  the  standard  state.  The  activity  now  possesses  a  peculiar  quality:  if  graphite  is 
present,  oq  =  1;  if  it  is  absent,  however,  it  is  not  true  that  is  zero,  nor  is  it  undefined. 
Consider  the  following: 


aco2  kco2  Vo2 


(A-3) 


When  graphite  is  present,  ap o2  depends  on  ao2  (or  vice  versa).  But  if  graphite  is  not  pres¬ 
ent  and  both  aco2  and  a02  are  known  (as  is  the  case  when  a  fuel  is  burned  in  a  large 
excess  of  oxygen)  the  activity  of  carbon  is  given  by 


aC  aC02  ^KC02  fl02  ) 

This  concept  has  been  used  throughout,  and  it  has  proved  a  useful  one. 


(A-4) 
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Appendix  B 

THE  NEWTON-RAPHSON  METHOD  FOR  SIMULTANEOUS 
NONLINEAR  EQUATIONS 


Suppose  the  variables  to  be  xj,  j  =  1  to  n,  given  n  nonlinear  equations  of  the  form 

fi(xJ,x2 - )  =  0,  i  =  1  to  n  (B-l) 

Let  xj  be  the  initial  approximation  for  each  variable,  from  which  can  now  be  computed 
the  set  of  fj(x°j  ,xi  ,  . .  .).  The  objective  is  to  improve  the  approximations  so  that  each 
fj  becomes  more  nearly  zero.  If  dxj  is  the  improvement  to  be  made  to  xj,  a  new  set  of 
approximations  is  generated,  xj'  =  xj  +  dxj  and  a  corresponding  new  set  of  fj.  Let  the 
change  in  fj  be  denoted  dfj  so  that  tj  =  +  dfj.  Since  fj  is  a  state  function  of  the  variables 

xj,  it  follows  from  the  principles  of  calculus  that 


(B-2) 


Hence  the  improved  set  of  equations  will  be 


If  j  D 

set  of  equations  has  the  matrix  form 


(B-3) 


the  vector  of  dxj,  Yj  the  vector  of  fj.  then  the 


D.. 

<5. 

= 

Y. 

U 

J 

l 

whence 

«J 

= 

-  Du  ‘ 1 

(B-4) 


Tlti-  latest  values  of  the  variables  are  used  to  compute  |  Djj  J  and  [  Y; 
evaluation  of  the  partial  derivatives  is  used.  I  6j 


;  a  numerical 


gives  the  next  approximation  from 


I  ||  '  ' 

which  a  new  set  of  Dy  and  Yj  are  found  and  the  process  is  repeated  to  any  desired 
precision.  1  1  J 
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In  this  case  each  fj  represented  the  error  in  the  material  balance  for  an  element.  Iteration 
was  continued  until  the  error  for  each  element  was  no  more  than  one  part  in  10^,  relative 
to  the  total  amount  of  the  element  present. 
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Appendix  C 

ISENTROPIC  PROCESSES  FOR  AN  IDEAL  GAS  OF  VARIABLE  COMPOSITION 


When  the  number  of  moles  of  gaseous  components  can  vary,  as  a  result  of  shifting 
chemical  equilibria,  the  ideal  gas  law  remains  P  =  nRT/V.  but  now  n  is  a  variable,  dependent 
on  P  and  T.  For  an  isen  tropic  process4 


dU  =  -PdV 


or 


CvdT  + 


dV 


-PdV 


(C-l) 


CV  = 


and  is  evaluated 


by  measuring  internal  energy  (relative  to  any  chosen  zero 


datum)  at  two  temperatures  which  are  very  close  together  and  ior  which  the  same  set  of 
condensed  phases  will  be  present  at  equilibrium.  The  partial  derivative  in  Equation  C-l  is 
found  from  the  thermodynamic  equation  of  state4 


(C-2) 


and,  for  variable  n. 


(C-3) 


When  Equations  C-2  and  C-3  are  substituted  into  Equation  C-l . 

cvdx  -  -  E  |„t  Tj|aj  Jdv 

that  is. 


In  T  =  - 


In  V 


(C-4) 
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An  “isentropic  coefficient”.  X,  is  now  defined  by  the  equation 


(C-S) 


so  that,  by  substitution  into  Equation  C-4, 

TV^  '  =  constant 

Cy,  as  defined  above,  and  X,  as  defined  by  Equation  C-5,  have  been  computed  for  each 
point  used  in  these  calculations.  Values  for  Cy  vary  widely,  being  especially  large  in 
regions  where  phase  changes  are  taking  place.  A  few  values  of  Cy  and  X  are  given  in  Table  8. 
(All  these  computed  values  have  been  retained  on  file  for  any  interested  person,  as  well 
as  values  of  the  average  molar  mass  and  the  number  of  moles  of  gaseous  products.) 
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